Inflationary cosmology with two-component fluid and thermodynamics 



Q"^ ' Edgard Gunzig [|, Alexei V. Nesteruk and Martin Stokley 

00 

Running head: Inflationary cosmology with two-component fluid and thermodynamics 

> 

m 

00 

o 
\o 
o, 

Q^ , Correspondence: 

o . 

qh, a V Nesteruk 

^ [ School of Computer Science and Mathematics 

^' Portsmouth University 

Mercantile House 
K/i I Hampshire Terrace 

'1-1 ■ Portsmouth POl 2EG 

^ ■ England 



TEL: + (01705) 843 108 
FAX: + (01705) 843 106 
EMAIL: alexei.nesteruk@port.ac.uk 



'^Universite Libre de Bruxelles, Service de Chimie-Physique, Campus Plaine CP 231, 1050 Bruxelles, 
Belgium 

School of Computer Science and Mathematics, Portsmouth University, POl 2EG England 



Abstract 



We present a simple and self-consistent cosmology with a phenomenological model of quantum 
creation of radiation and matter due to decay of the cosmological constant A. The decay 
drives a non-isentropic inflationary epoch, which exits smoothly to the radiation-dominated 
era, without reheating, and then evolves to the dust era. The initial vacuum for radiation and 
matter is a regular Minkowski vacuum. The created radiation and matter obeys standard 
thermodynamic laws, and the total entropy produced is consistent with the accepted value. 
This paper is an extension of the model with the decaying cosmological constant considered 
in Q . We compare our model with the quantum field theory approach to creation of particles 
in curved space. 



1 Introduction 

The aim of this paper is to extend the scenario of the evolution of the universe with smooth 
exit from inflation, and particle production at 

the expense of the decaying cosmological constant, developed in paper |j|]. In that paper 
thermodynamics and Einstein's equations led to an equation in which the Hubble rate H is 
determined by the particle number N. The model is completed by specifying the particle 
creation rate T = N /N, which led to a second-order evolution equation for H. The evolution 
equation for H then has a simple exact solution, in which a non-adiabatic inflationary era 
exits smoothly to the radiation era, without a reheating transition. For this solution, there 
were given exact expressions for the cosmic scale factor, energy density of radiation and 
vacuum, temperature and entropy. 

In the paper ^j we generalised the abovementioned results for the case of a scalar field 
ip interacting with radiation via the gravitational field, leading to cosmological evolution 
with smooth exit from inflation. Our particular task was to determine whether the theory 
formulated in terms of the scalar fleld (/? could lead to any new results in comparison with 
the previous case of a decaying cosmological constant. We concluded that the presence of 
the scalar field ip in this model did not change considerably the physical results which have 
been obtained in the paper ^. As result we argued that models with decaying cosmological 
constant A corresponding to a special case of the equation of state p^ = — P^p ■, describe 
adequately the smooth transition from inflation to radiation and give a reasonable prediction 
for the entropy of matter in the universe. 

In this paper we continue to argue along the lines of the paper ^ generalising its results 
for a two-component cosmological fluid, i.e. radiation and matter. The aim of the paper is 
to build a scenario of overall evolution of the universe with smooth exit from inflation to the 
radiation dominated stage and then its further evolution to the matter-dominated universe. 

The background of this paper is constituted by two ideas in physical cosmology. One 
of them is connected with the longstanding attempt to explain all matter in the universe 



as produced by quantum creation from vacuum. This has been studied via quantum field 
theory in curved spacetime (see for example [^, ^, |^, |6|, 0). Most cosmological models ex- 
hibit a singularity which presents difficulties for interpreting quantum effects, because all 
macroscopic parameters of created particles are infinite there. This leads to the problem of 
the initial vacuum (see discussion in [Q]). One attempt to overcome these problems is via 
incorporating the effect of particle creation into Einstein's field equations. For example, in 
the papers of the Brussels group Q , the quantum effect of particle creation is considered in 
the context of the thermodynamics of open systems, where it is interpreted as an additional 
negative pressure, which emerges from a re-interpretation of the energy-momentum tensor. 
This effect is irreversible in the sense that spacetime can produce matter, leading to growth 
of entropy, while the reverse process is thermodynamically forbidden. These results were 
generalized in a covariant form in Q. Our approach differs from that of |^, y] in that we 
do not modify the field equations. Instead, we associate the source of created particles as a 
decaying cosmological constant A. 



A number of decaying vacuum models has appeared in the literature (see [10, 11] and 
references cited there) . A review of the different phenomenological models of evolution with 
variable cosmological term can be found in paper p^ . 

Inflationary models with fixed cosmological constant and cold dark matter have been 
successful in accounting for the microwave background and large-scale structure observations. 



while also solving the age problem (see [13|). However, these models are challenged by the 
reduced upper limits on A arising from the Supernova Cosmology Project (see [0]), and 
also by the long-standing problem of reconciling the very large early-universe vacuum energy 
density with the very low late-universe limits jll[ |. 

One resolution of these problems is a decaying cosmological constant A which is treated as 
a dynamical parameter. This approach was typical for the quantum field theorists for many 
years (see for example ||l5[| ). Many potential sources of fluctuating vacuum energy have been 
identified which were to give rise to a negative energy density which grows with time, tending 
to cancel out any pre-existing positive cosmological term and drive the net value of A toward 
zero. Processes of this kind are among the most promising ways to resolve the longstanding 



cosmological 'constant' problem (see |16] for review). It is worth mentioning the recent paper 
of Parker |T^ indicating an attempt to revive the idea of the cosmological constant as a 
purely quantum effect associated with the renormalization of the general relativistic action. 

In ad hoc prescriptions, the functional form of A(t) or A(a) or A{H) (where o is the scale 
factor and H is the hubble rate) is effectively assumed a-priori (see the review ||l^ where all 
known forms for A are listed). Typically, the solutions arising from ad hoc prescriptions for 
A are rather complicated, and moreover, it is often difficult to provide a consistent simple 
interpretation of the features of particle creation, entropy and thermodynamics. 

In contrast to many other models, we propose a simple, exact and thermodynamically 
consistent cosmological history. The latter originates from a regular initial vacuum. Together 
with naturally defined asymptotic conditions for the number of created particles this leads 
to a simple expansion law and thermodynamic properties, and to a definite estimate for the 
total entropy in the universe. Since the exit from inflation to the radiation era is smooth, we 
avoid the problem of matching at the transition. A similar smooth evolution has been used 
in I, 111, 111, 1^. 

The choice of a as dynamical variable and the very simple form of H{a) that meets the 
physical conditions, lead to elegant expressions for all parameters describing the radiation 
and decaying vacuum, and also to a physically transparent interpretation of these results, 
including the estimate of entropy. 



We use units with SvrG, c and ks equal 1. 

2 Model 

Metrics and Matter 

We consider a spatially flat FRW universe we has the metric 



ds'^ = -df^ + a\t) dx^ + dy^ + dz^ , (1) 

containing a uniform two-component cosmological fluid, which is composed of non-interacting 
matter and radiation. The radiation has an energy density pr{t) and a pressure pr = Pr/'i 
while the non-relativistic matter has an energy density Pm{t) and pressure Pm = 0. (we 
assume that Pm ^ Pr-, nm *C n^ at high temperature 

r, so that we can neglect pm = n-mT with respect to pr). 

The energy momentum tensor of these components correspondingly is 

Decaying vacuum 

We also consider matter corresponding to the quantum vacuum energy, with energy momen- 
tum tensor 



SvrG 



g^lu 



(We shall use units such that SvrG = 1 = c henceforth) The conservation equations V'^(T^ -|- 
Tjf, + T%) = reduce to 

Pr+Pm + ^Hpr + ^H Pm = -A , (2) 

(which is a special form of the first law of thermodynamics), where H = a/ a is the Hubble 
rate. We can rewrite this equation, introducing an enthalpy of radiation and matter 

il — ilj' -}- ilrn 7 '^r ^ Pr i Pr — ^/ *J Pr •, '^m ^ Pm i Pm — Pm 

so that 

p + 3Hh = -A, (3) 

where p = Pr+Pm- The equations (S) and (H) show how energy is transferred from the vacuum 
to the radiation and matter densities. This energy transfer can therefore be understood as 
creation of quanta of radiation and matter from the vacuum. Indeed, employing the extended 
form of the first law of thermodynamics suggested in Q: 

d(pV)+pdV--d(nV) = 0, (4) 

n 

one can connect the evolution of p and p with the evolution of the total number of particles 
(both photons and massive particles) A^ = nV , where ^ is a comoving volume of the 
observable universe. Since (Q) is equivalent to 

p + 3Hh = h^, (5) 



comparing with (y) gives 

N A 

Therefore in order to create matter or radiation we need A < 0, i.e. A to decrease with time. 
It is clear from this formula that if A = const we have no particle production and the total 
number of particles is conserved. Since radiation and matter do not interact it is legitimate 
to assume in this case that radiation and matter evolve separately according to standard 
conservation law 

Pr + iHpr = 0, pm + SH p^ = 0, (7) 

which in conjunction with 

d{prV) + PrdV - —diUrV) = , d{p^V) + p„,dV - —diUmV) = , (8) 

tell us that the number of photons Nr = rirV and the number of massive particles N^ = nmV 
are conserved separately, so that N = N^ + N^ = constant. 

This result gives us a clear understanding that in the case of A = the initial vacuum 
for photons and massive particles (where Nr = Nm = 0) will be stable leading to no particle 
creation in the universe. By switching on the source A 7^ in the right-hand side of the equa- 
tions (^) and (^) we effectively switch on coupling of radiation and matter with gravitational 
field. In other words the evolving A(t) acts as an interaction of radiation and matter with 
the gravitational field leading, according to (^) to creation of photons and massive particles 
from vacuum. 

Field equations 

The field equations 



are 



3H' = pr + pm + A, (9) 

2H + 3H^ = --pr + A, (10) 

and if both are satisfied then the energy conservation equation (|2|) follows identically. 

Following |]l8|, we use a as a dynamic variable instead of t, and consider the Hubble rate 
as H = H{a) (in this case we cannot consider a = constant as a limiting case for a flat 
universe). Given H[a), we have two equations (|^) and (^ from which we can determine 
Pr{a),pm{a),A{a). 

Combining (^) and (^ one finds 

h{a) = -a[H\a)]', (11) 

where primes denote d/da. If we differentiate ([lO|)with respect to time and equate with (^), 
upon substitution for pm from (ll^) one finds an equation for p^. From this and (^) we also 
find an equation for pm 

Pria) = 3H'^{a)-R{a)+3A{a), p„(a) = i2(a) - 4A(a), (12) 



where 

R{a) = 3a[H^{a)]' + l2H^{a) (13) 

is a scalar curvature for metric (|l]). If we consider a = as an initial point of the evolution 
of the universe, it is reasonable to assume that the beginning of the evolution is a vacuum 
for radiation and matter, i.e. 

Pr(0)=0, Pm{0)=0, 

which gives an initial condition A(0) = 3Hq, Hq = H{0). Obvious physical conditions for the 
energy densities are 

Pr{a) > 0, pm{a) > 0. 

Using these conditions, equations (^) we can obtain the restriction for A which holds for 
any a: 

a [H\a)]'{a) < A(a) - 3H\a) <^a [H^{a)]' . (14) 

Assuming, that a > 0, and H{a) > (i.e. have expansion of the universe), one can 
conclude, that this inequality makes sense only if 

H'{a) <0, 

i.e. H Is a, decreasing function of a with the initial value Hq. Since A(0) = Si^o > ^^ have that 

A(a) -3H^{a) < 0, 

i.e. A(a) is a decreasing function. Therefore we can deduce from (14) the formula for A(a): 



A(a) = -f{a)/2a[H^{a)]' + 3H^{a) (15) 

where 7(a) is a continuous function with the range of values 

l<7{a)<2. 

On substituting (15) into ( |l^ ) one finds 

Pr{a) = ^^^1^ {-a[H\a)]'), pM = ^^^^ i-a[H\a)]'). (16) 

These equations satisfy the physical conditions placed on the energy densities and it follows 
that if 7 = 3/2, pm = 0, i.e. we have a universe containing pure radiation; and vice versa, 
and if 7 = 2, /9r = 0, i.e. we have only matter in the universe. 

Now we have the formulas for pr,Pm and A in terms of H{a) and 7(a). The next step 
hence is to make some reasonable assumptions about them. 

The form of ^(a) 

In order to make our model consistent with the fact that at present state of the universe 
radiation and matter evolve adiabatically, and matter dominates radiation, i.e. according to 
standard red shift law pr ~ a~^ and pm ~ a~^, we assume that at the late stage of evolution 

Pr{a) 3(2-7(0)) am 



p^{a) 2(27(0) - 3) 



where am is the value of the scale factor at time when the density of radiation and matter 
are equal {am ^ W^^ adecoupUng) ■ From this formula one can find 7(a) for a ^> am'- 



7(a) 



6(1 + a/a„ 
4 + 3o/a„ 



(17) 



An amazing feature of this formula is that it adequately describes a smooth evolution 
of 7 from a = to a = 00. Indeed 7(0) = 3/2, and lima^oo7(a) = 2. This result is in a 
coherence with the obvious physical expectation that the universe is dominated by radiation 
near the beginning of expansion, and the universe is dominated by matter at the late stage 
of its evolution. In other words the formula ([l7|) gives us a simple form for 7(a) for all values 
of a. 

Assuming now that the simplest form for 7(a) for any a is given by (0), and plugging 

[-a [H^{a)]'] , pm{a) = — Pr , 



T7| ) into (16) and (15) we find that 

3 



Pr{a) 
A{a) 



4 + 3a/a, 
3 {I + a/am) 
4 + 3a/ am 



[a[H\a)]']+3H\a) 



(18) 
(19) 



3 The law of evolution for H(a) 



The Form of H(a) in the vicinity a=0 

To write down the above equations in a form which depends purely on the scale factor a, 
we need to make predictions about possible form of H{a). To do this we bring into account 
two physical assumptions about the nature of the evolution of the universe in the vicinity 
a = and when a — > 00. It is by examining the behaviour of H(a) in the light of these two 
assumptions we can predict its possible form. 

Let us consider the case when a — > 0. In this case we can neglect a/am ^ 1 with respect 
to all 0(1)- quantities in the expressions (p^, (|l8|), and ( [191 ) (i.e. 7 « 3/2, corresponding to 
a purely radiation dominated universe as a ^ 0). 

Substituting (19) and (11) into (|6|) one obtains the equation for N when a — > 0: 



77 



5 




H' 




H"} 


— 


+ 




+ 




a 




H 




H' 



which integrates to 



-a'lH^y^'^' 



N = A\ 
where ^ is a constant. This expression for N can be rewritten as an equation for H: 

d 



da 



H\a) 



1 N^/^{a) 
2A ^^5 



(20) 



(21) 



We require that initially there is a standard Minkowski vacuum for radiation, so that 
N{a) -^ and n{a) ^ as o ^ 0. This implies the limiting behaviour 



N{a) ~a", a>3, as a ^ . 

It then follows from equations (p^) and (p2|) that PI 

H{a) — > constant , as a ^ . 



f22) 



(23) 



H(a) in the vicinity of exit from inflation 

This form of the asymptotic behavior of H near a = tells us that we have a type of 
expansion a la inflation. It cannot be an exact exponential inflation (with H = constant) 
for al a, because in this case we would have p^ = Pm = 0, A =constant, i.e. there would 
be an eternal stable false vacuum universe with no production of radiation and matter. 
Since the underlying motivation of our model is to obtain the observable figures for the 
energy-density and entropy of radiation and matter in the universe, one must assume that 
an initially inflationary universe will evolve into a present-state universe. Since inflation is 
usually understood as an expansion with acceleration, i.e. with a > Q, oi H + aH' > 0, and 
at present we have decelerated phase, there must be exit from inflation such that a = 0, or 
H(^ = —aeH'^,He = H{a(.). This makes possible to calculate the pr,Pm,-^ at a = Og in terms 
of H,: 



Pr{ae) 



A(ae 



6 



4 -|- Zae/ar, 

Qae/am 

4 -|- Soe/a^ 

6 (1 + Oe/fflr 

4 -|- 2>ae/ar, 



m 



Hi 



H. 



2 e 

dfri 

2 +3//2r 



— ^ Pr{ae) <- Pr{ae) 



2 



2 



(24) 
(25) 
(26) 



(We assumed here ag ~ 10 cm, am ~ 10^^ cm, and neglected ae/am ~ 10^^^ with respect 
to all quantities of the order 0(1)). Comparing ( ]2^ ) and (^) one can see that the creation of 
matter during the inflationary stage is damped with respect to the creation of photons with 
the amplitude 10~^^, the universe is therefore dominated by radiation at this stage with little 
matter. 

Assuming now that created radiation is a black-body radiation, one can estimate the total 
number of photons created at exit from inflation. Using (24) for a standard calculation of 
entropy (see [18|) one can show, that Nr{ae) ~ 10^*^. 



H{a) for large a 

In order to find an asymptotic for H{a) for large a we employ again the assumption about 



aa 



an adiabatic evolution of radiation and matter at present, i.e. pr 
Using the formula (18) one can obtain an equation for H^{a) for large a: 



and p„ 



(5a- 



d 
da 



H\a) 



Aam + 3a 
3aam. 



a 
7a 



a 

SOrr 



4:ar. 



+ 



It is occurring that the De Sitter stage of the evolution of the universe appears inevitably for any value 
of 3/2 < 7 < 2 taken as a — > 0, because, as follows from (uW and IRH), A dominates evolution. One can show 
formally, that taking 7 — constant in dly), one can find that H^{a) ~ ^(^"'-f'l/'y _|_ (j^ gg that H{a) ~ C as 
a ^ regardless of the value of 7. 



which integrates to 



H\a) 



a 

dCLry 



+ 



1 



(We disregard an arbitrary constant of integration because of an obvious condition for an 
open universe that /7 — > as a — > c»). 

Matching of two asymptotics for H{a) 

We have arrived to two asymptotic formulas for H'^{a): 

r'it „\ - . „„„„+ „„ „ , n „„A Tj'i 



H^ia) ~ const as a ^ 0, and H^(a) ~ — ^ -\ — ^ as a 

a^ a-^ 



oo 



Our intention now is to combine these two asymptotics in order to find a formula for H{a) 
which is valid for all a and describes a smooth transition from inflation to the adiabatic phase 
of expansion contingent upon independent evolution of radiation and matter. 

In accordance with the method of the papers [0, |2|, we try the formula for H{a) in the 
following form: 

Om. r 



H^{a) = C 



+ 



(27) 



where g • n = 4 and p ■ m = 3, as suggested by the asymptotic behaviour of H{a) examined 
above. 

We can easily fix the constants W and d^ by imposing the condition for the exit from 
inflation 2H^ = -OeiH^)' in d^. We obtain that 

bi = al, dP = al/2. 

(It is interesting to note, that the condition for exit contains only the products qn and pm.) 
The final choice of q,n,p,m we will base on the requirement to have at a ^ oo the 
corrections to both leading terms in ([2^), i.e. 1/a^ and 1/a^ to be 0(l/a^) and 0(l/a^), i.e. 
we take q = A^n = l,p = 3,m = 1. The constant C in ( ^7|) can be found from equating both 
sides at a = Oe, so that finally (compare with a similar result from [pO| and [pl|]) 



H\a) = Hi 



6ar 



3am + 4ae 



+ 



2at 



ai + a'^ am{al + 2a?) 



In terms of the dimensionless variable x = aja^ (|2q ) can be rewritten in the form: 



HHx) = HI 



1 + 4/3/i 



1 



1 + x^ 



+ 2/i 



1 



1 + 2x3 



(28) 



(29) 



where ^ = ae/am ~ 10^^^. 

This formula describes the overall expansion of the universe starting from inflation with 
H{x = Q) = \/2He. For x <1 the asymptotic for H^ has a form 



H' 



2m 



l + 2fi- Afix^ - x^ + 



1 + 4/3^ 
At exit H{ae) = H{x = 1) = Hf.- For x > 1 the asymptotic for m is 



H' 



2Hl 

1 + 4/3/x 



Ai ^ 1 


1 ^ 


1 


X'^ x^ 


"2x6 " 


"^ 



+ ... 



One can obtain an estimate for H^ at a = a^n (the time when the matter and radiation energy 
densities are equal), or x = ;U^^: 



4 Dynamics of radiation and matter 

Calculation of pr, Pm and A 

All further calculations will involve the expression 

d[H\x)] _ ml 



where 



such that 



-a[H\a)]' 
uj{x) = 



dx 



X 



+ 3^- 



1 + 4/3^ 



X 



oj{x) 



(1 + X4)2 "^(1 + 2x3) 



1 



cj(x = 0) = 0, u;(x = 1) = -(l + 4/3;u). 

with the asymptotics for the inflationary period, x < 1 

uj{x) « 3^x3 + x* - 12fix^ - 2x^ + . . . , 
and the asymptotic for expansion after inflation, x > 1 



3 /i 1 3 /i 2 



It is easy to obtain from here 



U![X = p 



rM • 



The expressions for pr, Pm and A have the form 



Pr{x) 

A(a) 



^^2 



1 + 4/3^ 



4 + 3x/i 



w(x), Pm{x) = Xppr{x), 



SHl 3(1 + x/i) 



1 + 4/3/x 4 + 3x;U 



uj{x) +3H'^{x). 



(30) 

(31) 



The leading terms in the asymptotics of these expressions for x < 1 have the form 

6H^ 



pr{x) 

A(a) 



1 + 4/3/u 
1 + 4/3^ 



[3/XX^ + (1 - 9/4^2 )x''] , Pm{x) = Xppr{x) , 

[1 + 2^ - 7px^ - (2 + 3/4^)x''] 



(32) 
(33) 



For the period from the exit from inflation to the time when the matter and radiation energy 
densities are equal i.e. for 1 < x < p~^ one finds 



Pr{x) 

A(x) 



24i72 



1 + 4/3^ 
6^2 



1 



4 + 3x^ 



3p_ , 1 

4x3 



+ 



Pm{x) =xppr{x) 



\ + xp 



l + (4/3)^ I l + (3/4)x;U 



3 /i 1 3 /i 

4 x3 x^ 4 x^ 



+ 



ii + J__ lii 
x3 X* 2 x^ 



(34) 
.(35) 



On substituting x = /i ^ one gets 

Pr{a = am) = Pm{a = am) = Prix = fJ.^^) ^ 6i/g/Lf^ = 6iJg 



15 15 

A(a = am) = A(x = p'^) w 6H^ y / = — 



(36) 



From the time when the energy densities were equal till the present and later (the matter 
dominated period) i.e. for 1 <^ px <^ x then we find the leading terms are: 



1 



Pr{x) PS QH^ — 



X^ 



Pm{x) « QHI^, 



^(^) - 'i^r 



These equations describe the standard red shift law pj. ^ a ^ 

and pm ~ a~^ and show the adequacy of our choice for the form of H^{a) above 

Particle Numbers and specific entropy 



Solving, by integration, the equations for N^ and Nm which follow from (S) 



Nl 



Pr a 



N' 

rn 

Nm. 



Pm _|_ 3 
Pm a. 



(37) 



one can get obvious expressions 



Nrja) 

Nriae) 



Prja) 

Pr{ae) 



3/4 r 



Nmja) 

Nm{ae) 



Pmja) 

Pm{ae) 



(38) 



One can easily determine Nr{ae) = Nr{x = 1) from Pr{ae) using e.g. the results of [| 
Finally 

16u;(x)T^/^ 



Nr{x) = Nr 



(39) 



_4 + 3/ij;_ 

Figure 1 shows the evolution of y{x) = Nr{x)/Nr e over the expansion of the universe The 
number of quanta of radiation increases to a maximum value which is given by its limit as 

X — > oo 

A^.oo = Nre2V2, (40) 

(we used p <^1). 

We can not estimate Nm{o-e) from pm{ae), because we do not know the mass of the 
particles constituting matter. Instead of postulating this mass, we can connect Nm{a-e) with 
Nr{ae) through the asymptotic condition for specific entropy s 



s = lim ^^ , , 
3;->oo Nm{x) 



10®. 



(41) 



Using Pm = {a/am)pr, one can easily prove that 



Nmjx) 

-'•^ m. p 



'Nrjx) 
Nr, 



4/3 



16uj{x) 
A + 3px 



x\ 



10 



Taking limit x ^ cxd in the last formula one finds 



Nrr 



47V„ 



which together with (|4^) and (^Tj) gives 



Nr, 



V2 



Nr,_. 



From the expressions (^) and ( pSD one can obtain the asymptotics for the number of 
particles: 



Nr(a) ~ a 



21/4 



N^{a) ~ a' 



as a ^ 



For the photon creation rate (|37| ) for x = a/ag < 1 one finds from ( |32|) (using // = 10 ^^ 
and neglecting it with respect to 0(1) terms) we find 

iV^ Af,. 4 ^ ^ V 3/n + X ■ 



In the limiting case x — > we have 



Then from (p7|) we find that 

r — 

and 

For the specific entropy 



T,(x ^ 0) = ^^e 



N:^{a) _4iV;(a) _ ^^^^ ^21// + 8x. 



Nm{a) 3 Nr{a) ^ 'V Sfi + x 

Tm{x ^ 0) = 7^i/e 



(42) 



s(a) 



iV^(a) 



A^m(a) 
one can easily obtain the equation for its evolution using (|4^): 



s 


'Nr{a) 
Nr{a) 


Nn^ia)' 


1 Nr{a) 


s 


Nm{a)_ 


3 iV^(a) 



1 



r^ < 0. 



The decrease of the s is connected with its definition and with the fact that N^ evolves faster 
to than Nr- Since F^ — > as a — > oo we have 



s — > as a — > oo 
One can find the formula for the evolution of s{x) 

16w(x)T-i/^ 



s[x) 



4 + 3/xx, 






B*V2. 
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5 Particle creation in decaying cosmology and quantum field 
theory in curved space 

In this final section we compare our phenomenological description of matter creation in an 
expanding universe with some known results from quantum field theory in curved space. 
Namely we compare the results for the N^ presented in a previous section with the number 
of particles Nq created from vacuum through quantum effects (see e.g. 0, |5|, Q) 

—^ = caityWit) or — -^ = c— — -V^, 

where R is the Ricci curvature given by the formula ([T^), and the constant "c" takes different 
values for different fields. This formula describes creation of massless non-conformal particles. 
We argued in the paper [^ that this formula can be integrated, so that the total number of 
particles created from the initial state at x = to some x can be presented in the form 
(x = a/oe): 

r x'^R^ix') , , 



or 



Jo H(x') "^ 

where Nq^, is the number of particles created to a = Ce {x = \). If we assume now that our 
species of particles can be modelled by quanta, let say, of a minimally coupled scalar field, one 



can apply ( p3D for the universe which evolution is described effectively by the formula (2£). 
In this case we tacitly assume that our phenomenological description of transfer of energy 
from decaying A to radiation and matter can be alternatively understood as quantum particle 
creation from vacuum in the metric (Q) with the expansion law (^). This assumption makes 
sense since, as we have seen above, the most of particles is created during the inflationary 
stage, where the contribution from massive particles can be neglected. 

Comparing two graphs for y{x) = Nr{x)/Nre (formula ( |39| ) and z{x) = Nq{x)/Nqe (fo- 
mula ( |4^ ) presented at the Fig. 1, we clearly see that the behavior of z{x) and y{x) is quite 
similar and the asymptotic values of z and y are of the same order: 

^(oo) = ^ = 1.113.. y{oo) = §^ = 2^2, 



Since Nqe is a free parameter in th e formula ( |43| ) one can choose it from a natural condition 
that 



so that 



-' 'qoo — -'^roo 



Nqe = 2.533.. Nr 



We can argue now that our phenomenological theory of production of radiation and matter 
via smooth evolution of the universe from inflation to the radiation-dominated stage and then 
to the dust stage, which is based on the form for the Hubble parameter (p9|), can be treated 
as a mimic of quantum particle creation from vacuum. This gives in a sense a microscopic 
justification for our phenomenological model. 
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Fig. 1 Graphs for 



z{x) 



Ngix) 



N, 



y[x) 



qe 



Nrix) 

~W7 



in dimensionless units where x = aja^- The number of particles created through quan- 
tum effects Nq is given by the formula (|4^); A^gg is its value at the time of exit from 
inflation. N^ and A^^e are related to photon creation due to decay of A and are given 
by the formula (p9|). It is clearly seen the that the behavior of z[x) and y{x) is quite 
similar and the asymptotic values of z and y ar of the same order: 



zyoo) — 



N, 



1.113. 



y{oo) 



Nr. 



qe 



Nr 



2^2, 



6 Conclusion 

We have considered a simple and thermodynamically consistent scenario encompassing the 
decay of the vacuum, the creation of radiation and matter, and a natural smooth transition 
from inflationary to radiation- and then matter-dominated expansion. In order to treat all 
matter in the universe as created from a Minkowski vacuum, we impose the physical condition 
that the number of particles is zero in the initial vacuum state a = 0. Together with requiring 
finite particle production in the observable universe, this constrains the form of the Hubble 
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rate (29), and gives an inflationary universe with smooth exit to the radiation era and then 
to the dust era. Using ( p9| ) we calculated then the energy density for radiation and matter, 
their partlicle numbers and specific entropy for the overall evolution of the universe. 

We argued that our phenomenological theory of production of radiation and matter via 
smooth evolution of the universe from inflation to the radiation-dominated stage and then 
to the dust stage, which is based on the form for the Hubble parameter (p9|), can be treated 
as a mimic of quantum particle creation from vacuum initial vacuum and give a microscopic 
justification for our model. 
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